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ABSTRACT 

Testing the standard Shakura-Sunyaev model of accretion is a challenging task because the central region of quasars where accretion 
takes place is unresolved with telescopes. The analysis of microlensing in gravitationally lensed quasars is one of the few techniques 
that can test this model, yielding to the measurement of the size and of temperature profile of the accretion disc. We present spectro¬ 
scopic observations of the gravitationally lensed broad absorption line quasar H1413-I-117, which reveal partial microlensing of the 
continuum emission that appears to originate from two separated regions: a microlensed region, corresponding the compact accretion 
disc; and a non-microlensed region, more extended and contributing to at least 30% of the total UV-continuum flux. Because this 
extended continuum is occulted by the broad absorption line clouds, it is not associated with the host galaxy, but rather with light 
scattered in the neighbourhood of the central engine. We measure the amplitude of microlensing of the compact continuum over the 
rest-frame wavelength range 1000-7000 A. Following a Bayesian scheme, we confront our measurements to microlensing simulations 
of an accretion disc with a temperature varying as T oc ^ jjjost likely source half-light radius of R 1/2 = 0.61 x lO'*" cm 

(i.e., 0.002 pc) at 0.18;rm, and a most-likely index of v = 0.4. The standard disc (v = 4/3) model is not ruled out by our data, and 
is found within the 95% confidence interval associated with our measurements. We demonstrate that, for H1413+117, the existence 
of an extended continuum in addition to the disc emission only has a small impact on the inferred disc parameters, and is unlikely to 
solve the tension between the microlensing source size and standard disc sizes, as previously reported in the literature. 
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1. Introduction 


The central black holes in quasars accrete the matter in their sur¬ 
roundings and emit light over the whole electromagnetic range. 
The standard model of continuum emission from an accretion 
disc devised by Shakura & Sunyae^ ( |1973| l predicts that the tem¬ 
perature of the accretion disc, far from its inner edge, varies with 
the distance to the black hole as T oc Under the assump¬ 

tion of local blackbody emission, this implies that the emitted 


flux increases with the frequency as in fy 


/ 1/3 


The emission 


from the continuum is sometimes altered in outflowing gas, as 
revealed by the presence of deep absorption troughs observed 
in 20-40% of the quasars. Broad absorption line (BAL) quasars 
constitute a class of objects, which show absorption in high ion¬ 
ization lines (such as C iv) with velocities reaching several tens 
of thousands kilometers per second. The structure and kinemat¬ 
ics of this wind is complex, possibly constituted by multiple 
components (|Lamy & Hutsemekers 2004[|Borguet & Hutsemek 


[ers ]M0| [O^wd et al.|2015i r Variability studies indicate that 

the wind may be composed of high velocity clouds located at a 
few tenths of parsecs from the disc ( Capellupo et al.|[2013| l in 
some objects, or at distances of several (hundreds of) parsecs in 
others ( |Arav et al.|2013| l. 


* Based on observations made with ESO Telescopes at the Paranal 
Observatory (Chile). ESO program ID: 386.B-0337 


H1413-I-117 is a famous BAL quasar at redshift Zs - 2.556. 
The light rays originating from this object are bended by a fore¬ 
ground lensing galaxy at a putative redshift ziens ~ 1.0 (|Magain| 


et al. 


1988[ Kneib et al.||1998|l, resulting in the observation of 


four images of this distant quasar. Each of these images is ob¬ 
served through a different portion of the lensing galaxy, such that 
the stars in the lens act as secondary gravitational microlenses 
that selectively magnify the source on scales of about 0.01 pc. 
This scale matches roughly the expected rest-frame ultraviolet 
size of the accretion disc predicted by the standard model. 


Since 1989, (rest-frame) UV spectra of the individual lensed 
images have been obtained almost every five years. Those data 
have revealed that the continuum in image D is about 0.7 mag 
brighter than what is measured in the broad and narrow emission 
lines. This brightening of image D is due to microlensing mag- 


niflcation of the continuum (Angonin et al. 

19901 Hutsemekers 

1993 

|0stensen et al.||19971|Chae et al.|2001 

IFlutsemekers et al. 

2010 

hereafter Paper I). In 2011, we obtained new spectra of 


H1413-I-117 covering the UV-optical rest-frame domain (1000- 
7000 A). In the present paper, we use these spectra to study the 
chromatic change of microlensing. The data also reveal a differ¬ 
ent amplitude of microlensing between the unabsorbed and the 
absorbed continuum fluxes, suggesting two sources of contin¬ 
uum: a compact region that is microlensed, and a more extended 
region that is not (or less) microlensed. In addition, we estimate 
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how this extended continuum modifies the size and temperature 
profile of the disc, as derived from chromatic microlensing. 

The structure of the paper is as follows. In Sect.j^ we present 
the data. Sectionj^explains how we use the microlensing signal 
from our spectra to unveil the existence of a non-microlensed 
continuum, and to measure the chromatic dependence of the 
amplitude of microlensing affecting the compact continuum. In 
Sect. after demonstrating the robustness of our detection of 
an extended continuum, we characterize its properties and dis¬ 
cuss its physical origin. The constraints we infer on the size and 
temperature profile of the accretion disc are presented in Sect.j^ 

2. Observations and data reduction 

Our principal data set are spectra obtained in 2011 with the 
FORS2 and SINFONI instruments mounted on the Very Large 
Telescope (PI; Hutsemekers, program ID;386.B-0337). These 
data are described hereafter. We complemented those data with 
archival optical spectra obtained on 1989-03-07 at the Canada- 
France-Hawai Telescope, and on 1993-06-23 and 2000-04-21-26 
with the Hubble Space Telescope (see Paper I for a complete de¬ 
scription). 

2.1. The SINFONI data 

H1413+117 has been observed on 2011-03-27 with the SIN¬ 
FONI H+K grating (R ~ 1500), dithering the object between 
four positions on the detector and two different position angles 
on the sky (0°and 90°). This resulted in eight exposures of the 
system, 600 s each. The sky spectrum was simultaneously ob¬ 
tained on each frame. The seeing was around 0.6". Reductions 
were performed using the Esorex-SINFONI pipeline (v2.2.9). 
The eight exposures were processed 2 by 2 to obtain four data 
cubes (wavelength-calibrated and sky-subtracted) with positive 
and negative images or spectra in each of them. The spectra of 
the individual components of the quadruple-lensed quasar were 
extracted by fitting four Gaussian profiles to the positive and neg¬ 
ative images of the quasar for each wavelength plane indepen¬ 
dently. Fitting Moffat profiles gave identical results. Telluric ab¬ 
sorption lines were corrected using the spectra of standard stars 
normalized by a blackbody. After a careful examination of the 
eight spectra of each quasar image, we selected the four best ex¬ 
posures that provided spectra in perfect agreement. 

2.2. The FORS2 data 

The FORS2 spectra were obtained by positioning a O'.'4-wide 
multi object spectroscopy (MOS) slit through the A - D images, 
and through the B - C images. The spectra of B and C are ig¬ 
nored in this analysis due to the small amount of microlensing in 
those images (i.e., about 0.25 mag). The grism 300V was used to 
cover the 3600-9000 A spectral range. Each observation of A-D 
was repeated at two epochs approximately separated by the time 
delay between the two images. The first data set was obtained 
on 2011-03-05 (hereafter epoch 1), and the second on 2011-03- 
31 (epoch 2). The spectra were obtained in the same range of 
airmasses and parallactic angle at epochs 1 and 2 (i.e., airmass 
1.404-1.281 at epoch 1 and 1.398-1.278 at epoch 2) such as they 
can be compared one-by-one directly without uncertain differen¬ 
tial extinction corrections. The observations were carried out in 
the EORS2 spectropolarimetric mode. Exposures of 660 s each 
were secured with the half-wave plate (HWP) rotated at the posi¬ 
tion angles 0°, 22.5°, 45°, and 67.5°. Eor each image and grism. 
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16 spectra (two orthogonal polarizations, four half-wave plate 
angles HWP, two epochs) were obtained and extracted by fitting 
two Moffat profiles along the slit. The analysis of the polariza¬ 
tion is presented in a separate paper (Hutsemekers et al., in prep). 
We corrected the spectra of epoch 2 for slit losses caused by see¬ 
ing differences to compare the spectra at epochs 1 and 2 . Eor that 
purpose, we used the spectra of three field stars observed simul¬ 
taneously to H1413+117. The comparison of the stellar spectra 
allowed us to verify the absence of chromatic changes at a level 
larger than 5% (from 3600 to 9000A), and to derive for each po¬ 
sition of the HWP a multiplicative correction factor. An average 
factor epoch2/epochl ^ 0.87 was found with a scatter of ~ 5%. 
A more significant corrective factor (epoch2/epochl ^ 0.65) was 
found for the HWP orientation of 67.5°due to the poorer seeing 
of the data obtained at epoch 2. Because of the larger uncer¬ 
tainties affecting that data set, we excluded it from the present 
analysis. 

3. Evidence for two continuum sources from the 
microiensing anaiysis 

Our analysis focuses on the microlensing in image D, which is 
compared to image A. The latter image did not show evidence for 
significant microlensing in the rest-frame UV wavelength range, 
over the last ten years (Paper I), and the current data confirm this 
finding. Although the intrinsic variability of H1413-I-117 is rel¬ 
atively small, ambiguity between differences caused by intrinsic 
variability and microlensing can be present. Eortunately, our new 
spectra have been obtained at two epochs separated by 25 days, 
matching the time delay Af = 23 + 4 days between images A- D 
( |Goicoechea & Shalyapin|2010[|Rathna Kumar et al.|2014| l. Un¬ 
less otherwise stated, when we present spectra of A and D in the 
(rest-frame) UV range, we use the spectra separated by the time 
delay, namely the spectrum of A obtained on 2011-03-05 and 
that of D obtained on 2011-03-31. 

3.1. Methodology 


We follow two complementary methods to characterize the mi¬ 
crolensing signal. On the one hand, the amplitude of microlens¬ 
ing as a function of wavelength is measured using the “MmD” 


decomposition technique (|Sluse et al.|2007||F 

utsemekers et al. 

2010| Sluse et al. 2012[ |Braibant et al. 2014 

', briefly summa- 
, measurement of 

rized hereafter (Sect. |3. 1. 1)1. On the other hanc 


microlensing in the Civ absorption trough, is realized by cal¬ 
culating the equivalent width (EW) of the absorption line over 
small velocity slices (Sect. [LL2| l. 

3.1.1. The MmD 

Assuming that the observed spectra E,- of the lensed image i is 
the sum of a spectrum Fm, which is only macrolensed, and of a 
spectrum Fm/j, which is both macro- and microlensed, it is pos¬ 
sible to extract components Fm and FMfi using pairs of observed 
spectra. Indeed, considering a line profile, we can write 

F\ - MFm + MpFMfi (1) 

Fi - Fm + Fmh , ( 2 ) 

where M - M 1 /M 2 is the macroamplification ratio between im¬ 
ages 1 and 2 and p the microamplification factor of image 1. 
We assume that image 2 is not affected by microlensing. We can 
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Fig. 1. Spectra of images A (black) and D (red) in the spectral range 1300-1700 A, covering the Siiv-Civ lines. The left and middle panels show 
the spectra as observed in 1989, 1993, 2000, and 2011, after normalization to unity in the continuum. The three shaded areas correspond to regions 
where the absorption changed significantly over the past 20 years, as discussed in Sect. |3.2[ The top and bottom right panels show the change of 
intrinsic absorption over time for images A and D. 


conveniently rewrite these equation to isolate the components of 
interest: 


Fm 

Fmii 


-A 

IF, 

A-M 

1 A 

M 

IF, 

A-M 

\M 



(3) 

(4) 


where A = Mfi. Up to a scaling factor, Fm only depends on A, 
while FMfi only depends on M. A is the scaling factor between 
the Fi and F 2 continua. It can be accurately determined as the 
value for which Fm{A) = 0 in the continuum. The microamplifi¬ 
cation factor of the continuum is then estimated using fi - AjM. 

For the MmD of images A and D presented hereafter, we fix 
M - 0.40 + 0.02 prior to the decomposition (instead of the em¬ 
pirical determination performed in earlier papers). This value is 
consistent with values retrieved around Hfi + [Out] with SIN- 
FONl data (Paper M - 0.38 + 0.01), and with the mid-lR flux 
ratios obtained by Mac Leod et al. ( 2009 [ M - 0.40 ± 0.06). 


3.1.2. The equivalent width 

The EW, because it involves a normalization by the continuum, 
is a simple method to unveil differential microlensing between 
the continuum and another source of emission. By definition, 
the equivalent width of a line (in emission or absorption) is 


EW ^ 


-I>9) 


dA, 


(5) 


where 'F is the continuum intensity underlying the absorption or 
emission line and usually approximated by the continuum emis¬ 
sion on either side of the studied feature. The flux over the wave¬ 
length range of interest [Ai,A 2 ] is denoted F^. It is easy to see 
that the EW is only modified due to microlensing if F^ and (F 


are microlensed at different levels, i.e., if the microlensed fluxes 
are F' - and F' = the EW only changes if + iic- 


3.2. Microlensing around the C iv line 

Eirst, we focus on the region of the Civ/115491ine. Figure [T] 
shows the variations observed in images A and D between 1989 
and 2011. The visual inspection of the spectra, normalized in the 
spectral range 1445-1455 A, free of line emission, shows that the 
two spectra do not match, and, in particular, the flux in the emis¬ 
sion lines (e.g., the peaks of the Civ and Siiv-Oiv) is smaller 
in image D than in A. This has been observed previously and is 
caused by microlensing amplification of the continuum of image 
D. Indeed, because we have normalized the spectra in the con¬ 
tinuum region, the magnification of the continuum of D leads to 
an apparently less intense emission line in that image, because 
the line is actually unaffected by microlensing. 

The other important feature to notice in Fig. [T] is the differ¬ 
ence of absorption depth in A and D on the blue side of C iv. 
This is particularly noticeable in the range 1455-1505 A, delim¬ 
ited with three shaded areas, and where the absorption is less 
deep. On the one hand, we see a decrease of the depth of the 
absorption with time. This corresponds to a long-term intrinsic 
decrease of the absorption. On the other hand, we observe a rel¬ 
ative change of the absorption in A and D over time: prior to 
2000, the absorption is deeper in image D than in A; in 2000, it 
is similar in A and D; and finally in 2011, it is deeper in image 
A. Because of this inversion of the relative absorption between A 
and D in 2011, at that epoch it is impossible to perform a MmD 
with a zero continuum everywhere in Fm, as shown in Fig.[^ In¬ 
deed, when we enforce the continuum to be zero in Fm, setting 
fi - 1.75, we violate the positivity constraint on Fm at the level 
of the absorber in the blue side of the C iv line (see light gray line 
in Fig. ig. In order to get Fm = 0 in that range, it is necessary 
to allow for a larger amount of microlensing of the continuum: 
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Fig. 2. MmD performed on images A and D of H1413+117 in 2011, us¬ 
ing M = 0.4 (Sect. |3.3^ . The dashed blue line is the flux of image A. The 
thick black line is the fraction of the flux Fm unaffected by microlensing 
for p = 2.0. The red line is the corresponding microlensed fraction of 
the flux F (shifted up to ease legibility), the dotted horizontal line is 
the line corresponding to Fm^ = 0. The gray lines show the decompo¬ 
sition for two other values of p (see Sect. |3.2^ . The three shaded areas 
correspond to the three regions A1 - A2 - A3 as discussed in the text. 


p - 2.0 {Fm shown with a thick black line in Fig.|^. This new 
value of p yields Fm > 0 in the unabsorbed continuum flux. At 
the same time, a large fraction of the continuum flux is still visi¬ 
ble in FMfi (i-C-, the continuum emission above the dotted line in 
Fig.|g. Therefore, the MmD suggests that there are two compo¬ 
nents producing the continuum, one originating from a region, 
which is sufficiently compact to be microlensed, visible in Fmii, 
and another one that we detect in Fm, and which comes from a 
nonmicrolensed, and, hence, more extended region. We empha¬ 
size that this detection of an extended continuum emission with 
the MmD also implies that the compact and the extended con- 
tinua are differently absorbed. More specihcally, the extended 
continuum has to be more absorbed than the compact continuum 
(in the range 1455-1505 A). As explained in Appendix]^ if both 
continua were similarly absorbed by the outflowing gas, or if the 
most compact region was more absorbed than the extended re¬ 
gion, the flux from the latter would remain undetected with the 
MmD when imposing Fm = 0 at the wavelength of the C iv ab¬ 
sorber (i.e., A ~ 1500 A). On the other hand, we may not exclude 
larger magnification of the compact continuum. Fig.j^shows the 
decomposition for p — 2.5 (top gray line). In that case, the inten¬ 
sity of the extended continuum emission equals the intensity of 
the compact (microlensed) continuum. 

To better characterize this differential signal in A and D, we 
show in Fig. the time variation of the EW of the absorption 
in three different wavelength ranges (Fig. and [^; the range 
A1 = [1496, 1505] A ([-10274, -8532] km/s), is a flat absorption 
plateau, which barely changes profile over time, but whose depth 
is signihcantly decreasing; the range A2 = [1482, 1496] A ([- 
12984, -10274] km/s) covers the highest velocity part of the ab¬ 
sorption as observed in 2011; and Anally A3 = [1455, 1482] A([- 
18209, -12984] km/s) is the range where absorption was detected 
prior to 2011. The variation of EW, shown in Eig.|^ allows us to 
assess the signiflcance of the relative absorption differences de¬ 
tected in Eig. [T] The EW declines over the 22 years time span of 
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year 


Fig. 3. Time variation of the equivalent width of the absorption, in the 
wavelength ranges A1-A2-A3 (see legend), for image A (solid) and D 
(dashed). 


the observations, in agreement with the intrinsic decrease of the 
absorption depth seen in Eig.[T] In addition, while EWa < EWq 
before 2000, the reverse behavior is detected in 2011. The ob¬ 
served differences of EW in A1 and A2 are highly signiflcant. 
We demonstrate in Appendix that it is only possible to ob¬ 
serve EWa > EWd in presence of two sources of continuum 
light, which are differently absorbed and microlensed. The dif¬ 
ference of EW observed between A and D at other epochs can 
also be understood under this hypothesis, but in that case the in¬ 
terpretation is not unique, as EWa < EWd can also be explained 
by superimposed broad line emission at the wavelength of the 
absorber. 

To summarize, the MmD decomposition and the EW varia¬ 
tion reveal the existence of two sources of continuum, a compact 
continuum that is microlensed, and a continuum that is more ex¬ 
tended, and hence not microlensed. The MmD reveals the ex¬ 
tended continuum in a visual and intuitive way, while the EW 
allows us to assess the signiflcance of the observed signal. 


3.3. Microlensing around other lines: chromatic microlensing 

The spectra obtained in 2011 cover a large fraction of the wave¬ 
length range [1100, 7000] A (rest frame). Using the MmD tech¬ 
nique (Sect. |3.1.1| l, we can derive the amplitude of microlens¬ 
ing as a function of wavelength. Because of the lack of apparent 
chromaticity of M empirically measured in past optical and NIR 
spectra (Paper I), we may assume that there is no signiflcant dif¬ 
ferential extinction between A and 4D The measurement of the 
amplitude of microlensing, in the vicinity of the main emission 
lines, is repeated for each individual spectra. The standard de¬ 
viation between measurements (from three spectra for EORS2 
data and four spectra for SINEONI) is quadratically added to the 
uncertainty in the relative flux calibration between epochs (UV 
data), and used to calculate the standard error on the mean. The 


Munoz et al. 


(2011 1 suggests possible extinction in A, but the time 


variability of the various flux ratios does not allow us to assess the 
uniqueness of this interpretation. As pointed out by |Munoz et al.| ( (2011[ >, 
extinction estimates are very uncertain in H1413-)-l 17. However, extinc¬ 
tion in A is smaller than in B, hence reducing the observed flux of A by 
less than 10% (cf. Paper 1 for extinction in B). This has a negligible 
impact on our source size estimates (Sect.|^. 
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results of this procedure are summarized in Table [T] The MmD 
of the main emission lines is presented in Fig. As explained 
above, the measurements in the rest-frame UV account for the 
time delay. Although intrinsic variability of the absorption lines 
is noticeable between the two epochs, it mostly affects the value 
of fi derived from C iv. In that case, measurements of yu on spec¬ 
tra obtained at a single epoch nevertheless agree at the 3cr level 
with those presented in Table We may then be conhdent that 
the values of yu derived at a single epoch from the Hj6 and Ha 
lines are accurate, in particular, since they are identical within 
the errors to those derived in Paper I. 

The extended continuum detected around C iv (Sect. |3.2| l is 
also tentatively detected around Ly a -N v. In that case, the ev¬ 
idence for an extended continuum flux is set by the absorbers 
in the range [1170-1190] A (rest frame), which have been identi¬ 
fied as Hi absorption at z ~ 2.4 - 2.45 (line number 48 to 51 in 
Monier et al.|1998) l. They are therefore likely to be high velocity 
absorptions associated with H1413-I-117. Note that in the MmD 
of the Ly aline shown in Fig.|^ the positivity constraint in Fm is 
(marginally) violated around A ~ 1165 A. This is acceptable be¬ 
cause this absorption feature is due to intervening C iv (line num¬ 
ber 45, 46 in |Monier et al.|1998') l, and therefore may be intrinsi¬ 
cally different in A and D. On the other hand, the absorption in 
the blue side of O iv-Si iv sets almost no constraint on the pres¬ 
ence of an extended continuum in that wavelength range. This 
is because the well-identified high velocity absorption A1 - A2, 
visible in C iv, is shallower and less clearly detected, in the O iv- 
Siiv blend. 


Table 1. Amplification factors for the (D,A) pair, with fixed M. 


Lines 

t(A) 

A 

M 

P 

Ly a-i-N v 

1216 

0.85+0.01 

0.40+0.02 

2.12+0.11 

Si iv-hO IV 

1402 

0.70+0.01 

0.40+0.02 

1.74+0.09 

Civ 

1549 

0.80+0.01 

0.40+0.02 

2.00+0.10 

Cm] 

1909 

0.68+0.01 

0.40+0.02 

1.74+0.09 

H/ 

4861 

0.62+0.01 

0.40+0.02 

1.55+0.08 

Ha 

6563 

0.62+0.01 

0.40+0.02 

1.54+0.08 


The measurements for UV lines are performed on pairs of spectra sep¬ 
arated by the time delay (i.e., 2011-03-05 for A and 2011-03-31 for D). 
The NIR measurements are performed on single epoch spectra obtained 
on 2011-03-27. For Lya-l-Nv(resp. Civ), the values are obtained for 
Fm = 0 in the range [1170-1190] A (resp. [1496-1505]A). ForSiiv+Oiv 
and C III] , the value of yu is likely underestimated because of the un¬ 
known contribution of the extended continuum at those wavelengths. 

3.5 pc on the radial distance of the cloud to the black hole. At this 
distance, the separation between A and D corresponds to a lin¬ 
ear beam separation of 3.2 x 10'^ cm, (10^^ pc) hence about ten 
times smaller than the gravitational radius of the source, ensur¬ 
ing that the two lines of sight cross the same absorber. Therefore, 
the possibility that the light from image A and D intercept differ¬ 
ent clouds (or regions of a cloud with different opacities) is very 
unlikely. 


3.4. Differences of light paths between images 

Up to now, we have shown that the continuum flux from 
H1413 h- 1 17 originates from two regions, a compact one, which 
is small enough to be microlensed, and a more extended one, 
which is not microlensed. An alternative interpretation could be 
that significant inhomogeneities of the outflowing absorbing ma¬ 
terial lead to different absorption for each of the lensed images. 
This can be tested by estimating the linear comoving separation 
between the two different lensed images A and D as a function of 
the distance to the central black hole. The comoving linear beam 
separation for z > ziens ( Smette et al.|1992 Chartas et al.|2007] l 
is given by 


S(Za)^0 


DoiDs. 




(6) 


ICapellupo et al.|2013l |Wildy & Goad|2015[|Grier et al.|2015| l. 


where 0 is the separation between the lensed images, D is the 
angular diameter distance, and the subscripts o, I, s, a refer to 
the observer, lens, source, and absorber, respectively. 

In order to see if the comoving distance between light paths 
is larger than a cloud size, we need to assume a reasonable cloud 
distance to the central engine. Important constraints on the dis¬ 
tance and geometry of absorbing troughs have been set by vari¬ 
ability studies of BALs. Two scenarios are generally considered. 
One explaining the variability as caused by rapid crossing of an 
absorbing cloud, and one as a rapid change of ionization in the 
ionized material (e.g., Lundgren et al.|2007| Filiz Ak et al. 2012[ 


Recently, Capellupo et al. ( 2014| l carried out a detailed study of 
the variability of the A1-A2 absorbers in H1413 h- 1 17. Based on 
the comparison of the C iv and P v absorptions, they showed that 
the C IV absorption is saturated such that its variability is most 
likely caused by rapid crossing of an absorbing cloud. Assum¬ 
ing a so-called crossing-disc model, they derive an upper limit of 


4. Characteristics of the extended continuum 

In this section, we first try to characterize some of the proper¬ 
ties of the extended continuum, in particular its wavelength de¬ 
pendence, and its possible time variability. We then discuss its 
origin. 


4. 1. Wavelength dependence 

The MmD yields the measurement of the relative fractions of 
compact (cc) and extended (ce) continua from the total flux. If 
we define / as the fraction of the total continuum flux (F = + 

Fee, which is extended (i.e., / = Fcel'F), we find based on the 
MmD around Civ and Lya that / = 0.3 (i.e., 30% of the total 
flux is nonmicrolensed). This value of / implies that FedFcc - 
//(I-/) = 0.42. Larger values of/may be acceptable (Sect.|3.2 
Fig.0. 

Although we do not observe a significant difference of the 
amount of extended continuum between the Ly a and C iv wave¬ 
lengths, we cannot exclude a small chromatic change. Since the 
Ly a and C iv absorbers used to derive p correspond to clouds at 
different velocities, the opacity of the absorber (i.e., Tcc and Tee, 
Appendix [A|) may be different in Ly a and C iv, leading to a co¬ 
incidentally similar measure of p around those two lines, but the 
effect is unlikely to be large. At wavelengths redder than C iv, 
the situation is more complex as the lack of intrinsic absorp¬ 
tion in that range precludes the detection of the extended con¬ 
tinuum. At the wavelength of the Balmer lines in the NIR, we 
measure p ~ 1.5, both in 2005 and 2011. This constancy of the 
amplitude of microlensing suggests that the contribution from a 
nonmicrolensed continuum is negligible in the redder part of the 
spectrum (i.e. ,A > 4800 A) because the extended continuum is 
stronger in 2011 than in 2005 (Sect. [4!^ . Consequently, we con¬ 
sider in the following analysis that there is no contribution of the 
extended continuum at those wavelengths. 
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Fig. 4. MmD for Ly a +N v (top left), C iii] (top right), H/3 (bottom left), H a (bottom right) emissions. The values of A = M x p derived for each 
line are based on the individual frames, while the MmD is shown for the sum of the individual spectra. 


4.2. Temporal variation 


Our detection of an extended continuum in 2011 has been pos¬ 
sible owing to the microlensing of image D and to the differen¬ 
tial absorption of the compact and extended continua in 2011, 
the extended continuum being more absorbed than the compact 
continuum. We may wonder whether the extended continuum 
was present at earlier epochs, or simply could not be detected 
because of the time variation of the relative absorption, an in¬ 
sufficient signal-to-noise ratio in past data, and/or changes of 
the microlensing configuration. Actually, the MmD presented in 
Fig. 6 of Paper I shows that Fm < 0 in part of the A2 - A3 
absorption for the data obtained in 2000. Our Fig. also re¬ 
veals that EWa > EWd for A2 and A3 in 2000, confirming 
(at least from A2) a significant detection of the extended contin¬ 
uum at that epoch. Accounting for this absorption in the MmD 
(i.e., enforcing Em = 0 at the location of the absorber) leads to 
FcelFcc ~ 20% and to p ~ 2.0 + 0.12. Although the detection 
is more noisy, this suggests that the extended continuum was in¬ 
deed present in those data but with a smaller contribution. At 
earlier epochs, we also detect differences of EW between A and 
D but with EWd > EWa. Although this is not an unambigu¬ 
ous signature of an extended continuum since emission from the 
BLR can also explain this signal (see Appendix [B]l, it is plausi¬ 
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ble to attribute the origin of this difference of EW to the extended 
continuum. Contrary to the situation observed in 2000 and 2011, 
the extended continuum has to be less absorbed than the compact 
one (Appendix [B]). Because we do not know the optical depth of 
the gas covering the compact and extended continua, we cannot 
derive FcelFcc with the MmD or based on the EW prior to 2000. 

Another way to approach the issue of time variability of 
FcelFcc is to look at the variation of the effective amplitude 
of microlensing over time. When calculating p' determined by 
making Fm = 0 in the unabsorbed continuum adjacent to 
Civ (see Appendix we find it to be roughly constant in 
1993/2000/2005 (p' ^ 2.0 + 0.2, Paper I). In 2011, p' sig¬ 
nificantly drops to 1.75 + 0.07. Under the hypothesis that the 
amplitude of microlensing does not significantly change with 
time, this indicates that a decrease of p' results from an increase 
of FcelFcc (i.e., Eq |A.5[ ). When we only measure the ampli¬ 
tude of microlensing for the compact continuum (i.e., allowing 
Fm + 0, cf Appendix [A), we find p — 2.0 + 0.1 in 2011 (Ta- 
ble[T]), which is consistent with the amplitude of microlensing in 
1993/2000/2005. Moreover, the micromagnification p of the HjS 
continuum did not change between 2005 and 2011, also sup¬ 
porting the constancy of p in the last decade. This absence of 
variation of the microlensing signal agrees with the small rela- 
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live transverse velocity of the lens inH1413+117 ( |Mosquera &| 
|Kochanek|201 l| l. 

In summary, we find evidence for extended continuum in the 
past data, but the data also suggest that the fraction of extended 
over compact continuum varies with time. 


4.3. A physical model 


The extended continuum unveiled by our data may simply be 
genuine flux from the host galaxy, or scattered emission by the 
host, similar to that unveiled by Borguet et al. ( |2008| l in Type 1, 
and by [Zakamska et ah] ( |20061 l in Type 2 AGNs. This interpre¬ 
tation is however very unlikely. Indeed, we have shown that at 
least 30% of the total continuum emission is extended. However, 
H1413 h- 117 is a BAL quasar, and most of the light is occulted 
by the optically thick absorption in the blue wing of C iv, with 
no more than a few percent of the flux detected at those wave¬ 
lengths. Although this absorber is at least at the distance of the 
BLR (because it occults the C iv emission), it is unlikely to cover 
the whole host galaxy. 

Broad absorption line quasars like H1413-I-117 are known to 
be intrinsically polarized because of scattering of the continuum 
flux in an extended region possibly located at a radius compara¬ 
ble to the BLR radius or in a polar wind (Lamy & Hutsemekers 
2004 Schmidt & Hines ||l999| l. It is legitimate to postulate that 
the same region is at the origin of the polarization and the source 
of extended continuum, especially since the observed fraction 
of extended continuum is consistent with the polarization de¬ 
gree observed in H1413-I-117 (Hutsemekers et al., in prep.). The 
scattering could be caused by dust, or by electrons. If dust is 
responsible for the scattering, we expect a wavelength depen¬ 
dence of the extended continuum as T This is consistent with 
the apparent absence of extended continuum in H/3, but implies 
an extended continuum emission 2.7 times stronger at the level 
of Ly a than around C iv. Another possibility is scattering by 
electrons. This does not yield a chromatic dependence of the ex¬ 
tended continuum, in agreement with our observation in the blue. 
Under this scenario, it remains to be seen how we can explain 
the lack of extended continuum at the reddest wavelengths, e.g., 
around HyS and H a. A combination of both dust and electron 
scattering is also possible. 

Alternatively, to explain the UV peak in quasar spectral en¬ 
ergy distribution (and other inconsistencies between the standard 
disc model and the SED of quasars), Lawrence ( |2012| l has pro¬ 
posed the existence of a population of dense clouds, which partly 
reprocesses the intrinsic continuum from the accretion disc. A 
false continuum is produced in the UV coming from an extended 
source. Since these clouds mostly reflect the inner continuum 
through electron scattering, we expect the ratio of the two con- 
tinua to be roughly constant in the Ly a-C iv spectral range, in 
agreement with our observations. This model also predicts a de¬ 
crease of the fraction of extended continuum at reddest wave¬ 
length, compatible with the absence of significant extended con¬ 
tinuum around H fi. In addition, in that model, UV variability is 
interpreted as changes of the cloud covering factor. It can occur 
independently of variations of the intrinsic continuum and can 
explain the change of Fed Fee between 2000 to 2011 (Sect. [4~2j l. 
In this model, the pseudocontinuum flux is suspected to orig¬ 
inate from cold clouds located at about 30Rsch (Schwarchild 
radii) from the continuum, which corresponds to a distance of 
0.037pc = 1.16 X 10^® cm ~ 0.6 t/q (for logiMBHlMo) - 9.12, 
and (M) = 0.3 Mq in Eq.[7]l. Even if the clouds were located in 
the range 30 - lOORsch, the extended pseudocontinuum would 


be of the order of 1-2 770 , and would therefore be sensibly mi- 
crolensed. 

In conclusion, we can rule out the host galaxy as the source 
of the extended continuum. Instead, we postulate that the latter 
is scattered light by dust or electrons. Scattering by dust or elec¬ 
trons could take place in the same region as that at the origin of 
the polarization in BAL quasars. Dust scattering yields a chro¬ 
matic dependence, while electron scattering does not. Therefore, 
dust may explain the apparently small amount of extended con¬ 
tinuum at the wavelength of Hy 6 , but implies that we underesti¬ 
mate the contribution of the extended continuum at the level of 
Ly a by almost a factor 3. Variations of microlensing and of the 
absorption in front of Civ and Lya may possibly help to test 
this scenario. Simple electron scattering fails to explain the de¬ 
crease of extended continuum at redder wavelengths. The model 
of blurred reflection of the continuum proposed by |Lawrenc^ 
( | 2012 | l accounts for most of the observational properties of the 
extended continuum we detect, except that the cold clouds re¬ 
sponsible for the "continuum-reflection" are expected to be lo¬ 
cated too close to the black hole to remain unaffected by mi¬ 
crolensing. 


5. Constraints on the accretion disc parameters 


We now use our measurements of the chromatic microlensing 
(Sect. |3.3| l to derive the size of the compact continuum (which 
originates from the accretion disc) as a function of wavelength. 
Eollowing various other works (e.g., |Eigenbrod et al.|2008[|Bate 
et al.|2008) Anguita et al.|2008 1 , we test an accretion disc model 
with a temperature profile T oc R This kind of model leads 
to a dependence of the source size cr as a function of wavelength 
cr oc A.''. This model has two parameters, the source size cr (con¬ 
verted into the half-light radius R 1 / 2 ) and the power-law index 
V, for which we derive the probability density function (PDE), 
following the Bayesian scheme described hereafter (Sect. HI}- 


5.1. Methodology 

To turn the chromatic variation of microlensing into a quantita¬ 
tive information on the disc, we need to compare this signal to 
microlensing simulations calculated for different sizes and tem¬ 
perature profiles of the accretion disc, as well as for different 
fractions of compact objects in the lensing galaxy. The data to 
which the microlensing simulations are compared are the mea¬ 
surements presented in Table excluding p measured from C iii] 
and Si iv due to the uncertainty in the fraction of extended con¬ 
tinuum at those wavelengths. It is important to keep in mind that 
those values of p are derived assuming that the fraction of flux 
from the extended continuum equals the minimum value needed 
to reproduce our data (i.e. ,FeelFec - 0.42, Sect. |4.1| l. 

We follow a Bayesian scheme similar to that developed in 
[Bate et ^ ( |2008 1. Throughout the analysis, the sizes are ex¬ 
pressed in units of microlensing Einstein radius 770 , i.e.. 


m 


4G(M) Do, As 


Doi 


2.03 X 10*® 


cm 


<M) 

1 0.3 Mp, 


(7) 


where (M) is the average mass of the stars in the lensing galaxy. 
The different steps of our analysis can be summarized as follows: 


i) We use the inverse ray-shooting code of |Wambsganss| 
( 1998[ ) to generate microlensing magnification pattern of 
00x100 77 q ~0.657x 0.657 pc^. We fix the value of the lo¬ 
cal convergence and shear at the location of image D to 
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(a:, y) - (0.576,0.68), ( |Pooley et al. 2012| l, and calculate 
maps for ten different fractions /» of the total matter density 
in form of compact objects, logarithmically distributed in 
the range /, e [0.01,1]. Because of the small dependence of 
microlensing properties on the mass function of microlenses 
( [Lewis & Irwin|1996[|Kochanek|20()4| ), we assume that all 
the microlenses have the same mass. We then convolve the 
maps with the luminosity profile of the source. We model 
the disc as a two-dimensional Gaussian profile. The exact 
shape of the profile has little influence on the results ( jMor-j 
jtonson et al.j 2005 j l, provided that profiles with the same 
half-light radius R 1/2 are considered. We consider 51 loga¬ 
rithmically distributed half-light radii e [0.025,2.5] ?7o. 

ii) We generate a list of 2x10^ random positions to be extracted 
from the maps, excluding 500 pixels ( 5770 ) at the border of 
the map. 

iii) Based on those positions, we extract the magnification for 
each source size, and construct a chromatic curve represent¬ 
ing the magnification as a function of the source size. 

iv) The wavelength of observations A, and size of the source are 
related through the equation 


( 8 ) 


o- = 0-01 — 
(do 


Our size estimates can be compared to expectation from the¬ 
ory. The radius at which the disc temperature equals the pho¬ 
ton energy, kT — hclA^f^^t is given by the following expression 
( [Poindexter & Kochanek|2010HBlackburne et al.|201 Ij ): 


/?i = 9.7xlO‘G — 

pm 


4/3 


(^r( 




1/3 


cm, 


(9) 


where rj is the accretion efficiency (assumed to be 77 = 0 . 1 ), Mbh 
the black hole mass, and L^dd is the Eddington ratio. An alter¬ 
native expression for the size which depends only on the magni- 
flcatio n corrected luminosity can also be derive d (Mor gan et al. 


2010 or based on Teff and M in Davis and Laor 2011 and Laor 


and Davis [ 2014 ) 1 : 


R, = 3.86 X 10* 


15 / ^opt,45 
COS(7) 


1/2 


(-) 

\pml 


4/3 


cm. 


( 10 ) 


where ctq is the disc size (in practice, the width of a Gaus¬ 
sian profile) at the reference wavelength do- This relation 
enables us to transform the dependence of the amplitude of 
microlensing with the size of the source into a wavelength 
dependence of the microlensing signal. We compare each 
simulated microlensing light curve to the data, and calcu¬ 
late a ^ I, (yUsim -jUobs)^/cr^ ^,j,^, where the sum is over 
all the observed wavelengths, and cr^,obs is the uncertainty 
in the estimated microlensing magnification at each wave¬ 
length. When necessary, we perform linear interpolation of 
the simulated microlensing signal to derive the amplitude of 
microlensing at the wavelength of the observations. 

v) We calculate the likelihood of each pair of parameter (cr, v), 
via summing the likelihood of each light curve estimated as 
T)(-Tobs|cr, V, /») = exp(-x^/2). This enables us to derive, for 
each fraction of compact matter /,, a two-dimensional map 
of the likelihood of the source parameters. 

vi) We perform a weighted sum of the likelihoods associated 
with each /», using the dark matter probability distribution 
derived by Pooley et al. ( 2012| l, based on the study of the X- 
ray microlensing in the system. We then convert the likeli¬ 
hood into probability distribution using the Bayes theorem. 


5.2. Results 

The 2-D probability density function d^PI(dRi/ 2 dv) resulting 
from our analysis is shown in Fig. This PDF is marginalized 
over the fraction of compact objects /» in the lensing galaxy, 
weighted by the dPjdf, estimated by [Pooley et al.[ ( |2012| l. Our 
most likely values of the disc size and power-law index of 
the temperature profiles are Ri/ 2 {Ao - 0.18/2777) = 0.3 770 = 
6.1x10*5 V<M) /0.3 Mo cm = 0.00197 '^{M) 10.3 Mq pc, and 
V = 0.4. However, those values significantly depend on our prior 
on the fraction of compact objects. This is illustrated in Fig. 
where we show the PDF for /?i /2 at a fixed value of v = 4/3 (i.e., 
the standard accretion disc value), and for various amounts of 
the compact matter fraction /». Those PDFs are compared to the 
case where we use the prior on /, from Pooley et al. ( [2012 i.e., 
as in Fig.|^. 


where Lopt ,45 is the luminosity at 4861A in units of 
10^5 ergs ' cm“^, and i is the inclination of the disc. Figure 
shows that the theory size (with Lopt ,45 = 5.22; Sluse et al.|2012 1 
is somewhat larger than the microlensing size (using (M) = 
0.3 Mq in Eq. |^. The presence of an extended continuum only 
has a small impact on this result. If we repeat our analysis using 
microlensing factors derived by forcing = 0 in the contin¬ 
uum (e.g., p' — 1.75 in Fig. [T] see also Appendix [A), we then 
And the PDF represented by the dashed contours in Fig.|^ which 
peaks at (/?i/ 2 , v) = (0.37 770 ,0.2). 

Because we cannot exclude that the amplitude of microlens¬ 
ing in the compact continuum is larger than the value used above 
(for Feel Fee ~ 0.4), we experimented with two alternative hy¬ 
potheses for the chromatic behavior of the extended continuum 
(Sect. [4T| l. On the one hand, we assumed that in the blue (i.e., 
up to Civ), Feel Fee ~ yielding an increase of the amplitude 
of microlensing in the blue, namely p - 2.5 ±0.12 for Civ and 
p - 2.56+0.13 for Ly a. On the other hand, we considered a con¬ 
stant fraction of extended continuum observed from C iv to H a 
(this yields an increase of the amplitude of microlensing p ~ 2 
around H/? and Ha). The most likely sizes derived for those al¬ 
ternative signals fell within the Icr solid contours of Fig. We 
can then conclude that our estimates of R 1/2 and v are robust 
within a factor 2 with respect to Fed Fee, and that the strongest 
source of error is the prior on /». 


5.3. Discussion 

In recent years, microlensing studies based on significant sam¬ 
ples of microlensed quasars found that the size of the accretion 
disc is larger by a factor 2 to 10 than the theory size derived from 
Eq. [9| ([Morgan et al.[[201(){ [Blackburne et al.[ 201 1[ [Jimenez-) 
Vicente et al. 2012[ and see Hall et al. 2014 for a summary). 
Our snapshot (i.e., single epoch) microlensing measurement of 
H1413+117 instead favors a most likely half-light radius smaller 
than the theoretical prediction, although the upper limit of the 
confidence interval extends to about two times the theory size 
(99.7% C.I.; Fig. ]^. There is therefore no tension between the 
disc size that we measure in H1413+117 and that measured in 
other systems. More interesting is the impact of the presence of 
an extended continuum amounting to 30% of the total flux. Fig¬ 
ure shows that the disc size shrinks by typically 20% when 


Morgan et al. ( [2010) l 


we account for the extended continuum, 
and [Dai et al. ( 2010) 1 reached similar conclusions, assuming that 
30% of the broadband flux in lensed quasars comes from regions 
of larger physical extent than the compact continuum. 
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Fig. 5. Probability distribution for the size of the disc 7?i/2, and for 
the power-law index v of the temperature profile. Solid contours cor¬ 
respond to the constraints obtained after correction of the microlensing 
amplitude for the presence of an extended continuum, and the dashed 
contours show the constraints if we do not account for the extended 
continuum. Contours are 68.3, 95.5, and 99.7% confidence level. The 
corresponding peaks of the prohability distribution are shown as a filled 
black star and red cross, respectively. 



Fig. 6. PDF of the source size for a standard accretion disc model 
(v = 4/3) and various fractions of compact objects /. (extreme values 
/,=!, shown in light gray, and /,=0.01, shown in dark gray, correspond 
to 100% and 1% of the total matter being compact, respectively). The 
thick blue line is the result obtained assuming the same prior on /,, as 
in Fig. 1^ namely the prior resulting from the analysis of |Pooley et al.| 
\2012). The vertical solid line is the theory size for a Shakura-Sunyaev 


model based on Eq. 10 for cos(i) = 0.5 (R 1/2 ~ 3.04 10'^ cm ~ 10 ^ pc). 


The other constraint provided by our data regards the power- 
law index of the disc temperature profile. We find a most likely 
index, v = 0.4, which is significantly smaller than that of the 
standard accretion disc, v = 4/3 (Fig.|^. Nevertheless, the lat¬ 
ter is compatible with our data within 2 cr, as generally found 


for measurements based on single epoch data (e.g., Bate et al. 
|2008[ Floyd et ^|2009| l. Our result also agrees with the find¬ 
ing of Blackburne et al. ( 201 l| l that microlensing measurements 
support a steeper temperature profile than in the standard model. 
Again, we should underline the rather small impact of the ex¬ 


tended continuum on those results. Conversely, the fraction of 
compact objects in the lens does play an important role on the 
parameter estimates. As suggested by Fig.[^ our data are more 
easily reproduced by a standard disc for higher values of /», al¬ 
though a large fraction of compact objects is not favored by X- 
ray microlensing studies ofH1413-i-117 ( Pooley et al.|2012) l. 

In summary, the extended continuum emission is not suffi¬ 
cient to explain discrepancies generally found between the stan¬ 
dard accretion disc model and the microlensing measurements, 
however, we emphasize the important role played by the prior 
on the dark matter fraction on size measurements. Dark matter 
fraction estimates result from analyses, which generally assume 
the size of the source. pimenez-Vicente et ^ ( |2015| l has recently 
demonstrated that this assumption has a non-negligible effect on 
the derived dark matter fraction at the location of the lensed im¬ 
ages. The common assumption of a very small source size (es¬ 
pecially in the X-rays) in dark matter fraction inferences yields a 
fraction of mass in compact form (/,) smaller than what is found 
by accounting for the extended nature of the source. If we pos¬ 
tulate that the prior on /» from Pooley et al. ( 2012| l is shifted 
to low values because of this source size effect, we then find 
that the most likely disc parameters derived from microlensing 
in H1413-I-117 are closer to those of the standard accretion disc 
model. 


6. Conclusions 


The technique of quasar microlensing, which consists of ana¬ 
lyzing the deformation of gravitationally lensed quasar spectra 
produced by stars in a lensing galaxy, is a powerful tool for 
studying the structure of the accretion disc since it constrains 
the size of the microlensed emission. This method generally 
shows that the size of the accretion disc is significantly larger 
than the theory size (see [Hall et al.|2014 for a good summary). 
In this work, we have analyzed new spectra (obtained in 2011) 
of the lensed images A and D of the cloverleaf H1413-1-117, a 
quadruple-imaged broad absorption line system located at red- 
shift z ~ 2.55. Because of the presence of absorption covering 
the source of continuum, and of the existence of microlensing 
in one of the lensed images, we have unveiled the presence of 
two sources of continuum: a compact continuum-emitting re¬ 
gion that is small enough to be microlensed, and a spatially sep¬ 
arated nonmicrolensed continuum-emitting region that accounts 
for about 30% of the total continuum in the UV rest-frame part 
of the spectra (namely 1000-3000A). Because this extended con¬ 
tinuum is detected in the UV, and is absorbed by the broad ab¬ 
sorption trough, we think that it may not come from the host 
galaxy, but from the innermost regions of the system (possibly 
as small as the quasar torus or the broad line region). This ex¬ 
tended continuum may be scattered light reprocessed either in 
the same scattering region as that at the origin of the polariza¬ 
tion in BAL quasars, or closer to the compact continuum source 
as in the model proposed by |Lawrenc^ ( |2012[ ). A reanalysis of 
archive spectra, shows that this extended continuum was already 
present in the past 22 years, but was less prominent compared to 
the compact continuum. 

Owing to our method of analysis of the spectra, we ef¬ 
fectively deblended the compact continuum from the nonmi¬ 
crolensed continuum. The amplitude of microlensing for the 
compact continuum, which we attribute to the accretion disc 
emission, decreases as the wavelength increases. We use this 
chromatic variation of the microlensing to derive the size and 
index of the temperature profile of the disc, assuming that the 
latter varies like T oc (where v = 4/3 corresponds to 
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the standard Shakura-Sunyaev value). For this purpose, we de¬ 
velop a Bayesian analysis scheme, which compares simulated 
microlensing signals for different source sizes and index v to the 
observed amplitude of microlensing. This analysis yields a most 
likely half-light radius of the source of R 1/2 - 0.61 x 10'® cm at 
0.18 pm, and a most likely index v = 0.4. The index is smaller 
than the prediction from the standard model of accretion, but re¬ 
mains compatible with the latter at the 2cr level. Fixing v = 4/3, 
we find a disc size that is about four times smaller than the theory 
size, but still statistically compatible. There is therefore no evi¬ 
dence for a stronger discrepancy between theory and microlens¬ 
ing disc size in this broad absorption line quasar compared to 
what is found in other systems. 

To test whether the extended continuum emission could ex¬ 
plain the tension between standard theory and microlensing re¬ 
sults, we repeated our measurement using an effective amplitude 
of microlensing, which is the amplitude of microlensing mea¬ 
sured including both the compact and the extended continuum. 
This is equivalent to the measurements performed for other sys¬ 
tems, where extended continuum emission is not detected (even 
when present). As expected, this procedure yields an increase 
of the microlensing size, and a decrease of the index v, but the 
change is not large enough to explain the large microlensing 
sizes that are commonly reported. We emphasize that another pa¬ 
rameter of the analysis, the fraction of matter in compact form in 
the lensing galaxy, may play a more important role in the source 
size measurement than the presence of an extended continuum. 
This is especially true for lensed systems where a small fraction 
of matter in compact form is assumed in simulations. 

Our detection of two spatially separated sources of contin¬ 
uum has been possible owing to differential absorption between 
the two sources of continuum, and to microlensing of the most 
compact continuum. The same method may be used in other 
broad absorption line systems, but its guarantee of success de¬ 
pends on the exact location of the absorbers with respect to 
the central engine. Innovative methods need to be developed 
to answer the many open questions about this extended con¬ 
tinuum such as: Is it present in other quasars, and visible for 
all quasar types? How does this continuum emission vary with 
wavelength? Is the extended continuum scattered light by the 
BAL wind, by the dust torus, or by cold clouds inner to the BLR? 
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Appendix A: MmD with an absorbed extended 
continuum 


We consider a line profile where absorption is present. We write 


Fi = M{Fee + F,ee-^-) + MfiF,,e-^-, (A.l) 

Fl = Fee + Fee + Fee (A.2) 


where Fee represents the flux from an extended source of con¬ 
tinuum, which is only macrolensed as the emission line flux Fee, 
Fee 5^ 0 is the flux from the source of continuum compact enough 
to be microlensed, and (Tee) is the optical depth of the ab¬ 
sorber in front of cc (ce). Using these relations in Eqs. 3 and 4 


with A = 

Mp, we obtain 



r A'-A 


II 

tin 

i'ee+i’ce^ + 1:4 ^ cc C 

A' - M 

(A.3) 


^ M-M' 


II 

t 

i' cc^ + ^ iF ee + i' ce^ ) ■ 

(A.4) 


In these equations, A' and M' are parameters tuned to satisfy 
constraints on Fm and F^/j- The constraint Fm/j = 0 at wave¬ 
lengths where Fee = 0 and Fee + 0 can only be satisfied 
with M' = M.\f M is correctly estimated, Fmi^ = Fee > 0 
independently of the existence of a nonmicrolensed extended 
continuum. 

In the absence of a nonmicrolensed extended continuum (i.e.. 
Fee - 0), making Fm = 0 at wavelengths where Fee - 0 and 
Tec — Tee — 0 (i-c., in the unabsorbed continuum) gives A' = A. 
In the presence of a nonmicrolensed extended continuum, mak¬ 
ing Fm = 0 in the unabsorbed continuum gives A' + A. In this 
case, Fm = Fee at wavelengths where the contin¬ 

uum is absorbed and Fee - 0, which can result in either positive 
or negative values for Fm whether Tee < Tec or Tee > Tec, respec¬ 
tively. Making Fm >0 allows the detection of the extended con¬ 
tinuum in the latter case (Tcc > Tcc), but A' = A is derived only 
if = 0 at some wavelengths where Fee - 0 and + 0. 
Without absorption (i.e.. Tee - Tee - 0), such a nonmicrolensed 
extended continuum would remain unnoticed. 

Note that /i' = A' jM, determined by making Fm = 0 in the 
unabsorbed continuum, can be seen as the average microampli¬ 
fication of the whole continuum source, while fi - AjM refers to 
the microamplification of the compact continuum source. They 
are related by 


- 1 Fee 


(A.5) 


Appendix B: Equivalent width in the absorption 
trough 

In order to understand how microlensing affects the EW, we can 
write Eq.[^ with more explicit expression of F^ and of the con¬ 
tinuum flux 'F. We assume that the latter is constant over the 
emission line and the absorption trough and is well approximated 
by the continuum flux measured in the vicinity of the line. We 
consider a broad emission line characterized by a flux F^^e, and 
a continuum originating from two regions, a compact region cc, 
compact enough to be microlensed, and a more extended region 
ce, which is not microlensed. We can specialize Eq.j^in the pres¬ 
ence of absorption. Eor a lensed image i, we obtain 


EW: = 


pAi 

Jai 


1 - 


e Fx,cc + e Fx,cc + e Fx,e 

E-i Fee + Fee 


dA, 


(BA) 

where /i, is the amplitude of micromagnification of the com¬ 
pact continuum cc, and Tee, Tee, and te are the optical depths 
of the gas absorbing the compact continuum, the extended con¬ 
tinuum, and the broad emission line, respectively. Contrary to 
Appendix [a] we make explicit the fact that the emission line can 
be absorbed (i.e.. Fee - F^^e)- 

Equation ED clearly shows that if F^ ee - F a,e - 0 (i.e., 
no extended continuum or emission line), and if the continuum 
does not vary over the range [Ai,A 2 ] (i.e., F^^ec ~ Fee), then 
the EW is independent of the amplitude of microlensing /r, . In 
the following, we compare the difference of equivalent widths 
AEW - EWd - EWa, between image D, characterized by mi¬ 
crolensing p, and image A, unaffected by microlensing, in the 
presence of flux from the emission line and/or from an extended 
continuum region. 


Appendix B. 1: Compact continuum+emission line 


We first consider the situation where there is only the compact 
source of continuum Fee the emission line over the wave¬ 
length range [di,/I 2 ]. We then have F^^ee = 0 in IB.ll Assuming 
that the continuum flux does not vary over the range [di, d 2 ] such 
that F^^ec - Fee, we can rewrite 


B.l 


EWi 


= (Pi Fee - Pi e Fee -e^’^ Fx,e) dA (B.2) 

- f (1 - dA - Fx^e) dA. (B.3) 


The difference of equivalent width, AEW - EWd - EWa, 
can be derived from the above equation specialized for the pair 
of images, i.e., 

u - 1 1 r^^- 

AEW^- -— {e-^'^Fx,E)dA. (B.4) 

P Fee Jxi 


It results from this equation that when p > 1, one may only 
observe EWd > EWa- 


Appendix B.2: Compact + extended continue 

We now consider that there is no flux from the emission line 
(Fx^e - 0) over the wavelength range [di, d 2 ], but two sources of 
continuum: a compact continuum microlensed by a factor /r, and 
an extended continuum Fee which is not microlensed. Assuming 
that the continuum flux is constant over the range [di,d 2 ] such 
that Fx,ec = Fee and Fx,ee = Fee, we obtain 

EWi = ^ (Pi Fee + Fee - Pi e’"" Fee - Fee) dA (B.5) 

= T i^. Fee (I- e-"“) + (e-^“ Fee)] dA, (B.6) 

where we defined F) = Pi Fee + Fee, the flux in the unabsorbed 
continuum. 

The difference of equivalent width AEW - EWd - EWa, is 
then 


A d^ce Fee " 1) j, ,r> -ta 

AEW = -——- (e “ - e “) dA. (B.7) 

J D J A JXi 

We clearly see that a difference of equivalent width is only 
detected if there is a difference of opacity between the absorbers 
covering the extended and the compact continua, i.e.. Tee ^ Tee. 
\f p> 1, we find AEW > 0 (EWd > EWa), and if Tec > Tee, and 
conversely AEW < 0 (EWa > EWd) if Tcc < Tee. 
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Appendix B.3: Compact + extended continue + emission iine 

The situation where the source of light under the absorption is 
the sum of a compact and an extended continuum and of an 
emission line, directly results from the two previous cases. If 
we define the flux of the unabsorbed continuum of image A and 
Fee + Fee, and To - p Fee + Fee, then we can write the 
following difference of equivalent width AEW - EWo - EWa 
between images D and A : 



From this equation, we see that when ju > 1 and Tec > Tee 
(i.e., when the compact continuum is more absorbed than the ex¬ 
tended continuum), KEW > 0. Conversely, tsEW < 0 can only 
be obtained if Tec < Tee, i-e., when the compact continuum is less 
absorbed than the extended continuum. Positive value of KEW 
may also occur depending on the absorption rate of the differ¬ 
ent components, and on the relative contribution of the extended 
continuum and of the emission line. 
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